Non-Markovian Persistence and Nonequilibrium Critical Dynamics 
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The persistence exponent 9 for the global order parameter, M(t), of a system quenched from the 
disordered phase to its critical point describes the probability, p(t) ~ t~ e , that M(t) does not change 
sign in the time interval t following the quench. We calculate 6 to 0(e 2 ) for model A of critical 
dynamics (and to order e for model C) and show that at this order M(t) is a non-Markov process. 
Consequently, 6 is a new exponent. The calculation is performed by expanding around a Markov 
process, using a simplified version of the perturbation theory recently introduced by Majumdar and 
Sire [Phys. Rev. Lett. 77, 1420 (1996)]. 
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The 'persistence exponent', 8, which characterizes the 
decay of the probability that a stochastic variable exceeds 
a threshold value (typically its mean value) throughout 
a time interval, has attracted a great deal of recent inter- 
est |l|dll) . One of the most surprising properties of this 
exponent is that its value is highly non-trivial even in 
simple systems. For example, 8 is irrational for the q > 2 
Potts model in one dimension p| (where the fraction of 
spins that have not changed their state in the time t after 
a quench to T = decays as t~ e ) and is apparently not 
a simple fraction for the diffusion equation |9],|l0| (where 
the fraction of space where the diffusion field has always 
exceeded its mean decays as t~ e ). 

A recent study of non-equilibrium model A critical 
dynamics, where a system coarsens at its critical point 
starting from a disordered initial condition, looked at the 
probability P(ti,t2) that the global magnetization does 
not change sign during the interval t\ < t < t 2 [Oj . 
The persistence exponent for this system is defined by 
P(ti,t2) ~ (ti/t2) e in the limit £ 2 Ai —> 00. Explicit re- 
sults were obtained for the ID Ising model, the n — > 00 
limit of the 0(n) model, and to order e = 4 — d near di- 
mension d = 4. For these systems it was found that the 
value of 8 was related to the dynamic critical exponent z, 
the static critical exponent 77, and 'nonequilibrium' expo- 
nent A (which describes the decay of the autocorrelation 
with the initial condition, (</>(x, £)</>(x, 0)) ~ £~ A / Z ) by 
the scaling relation 8z = A — d + 1 — 77/ 2. This relation 
may be derived from the assumption that the dynamics 
is Markovian, which is indeed the case for all of the cases 
considered in that paper. 

From a consideration of the structure of the diagrams 
which appear at order e 2 (and higher order), however, it 
was argued that the dynamics of the global order parame- 
ter should not be Markovian to all orders, implying that 
the exponent 8 does not obey exactly that 'Markovian 
scaling relation' Q|. This means that 8 is a new ex- 
ponent. Monte-Carlo simulations in 2 dimensions indeed 
suggest weak violation of the Markov scaling relation pij . 



In this Rapid Communication we present an explicit 
calculation of the non-Markovian properties of the global 
order parameter. The nonequilibrium magnetization- 
magnetization correlation function is calculated to order 
e , and this is then used to calculate 8 to the same or- 
der, using a perturbative method proposed by Majumdar 
and Sire (MS) ||, valid in the vicinity of a Markov pro- 
cess. The Markov scaling relation is shown explicitly to 
be violated at order e 2 , so 8 is indeed a new independent 
exponent. 

Before discussing the calculation of 8, however, we pro- 
vide first a simpler, and more transparent, formulation of 
the perturbation theory than that given in MS. In partic- 
ular the final result, Eq. (pi3), does not appear explicitly 
in MS (12). 

Let y(i) be a Gaussian stochastic process with zero 
mean, whose correlation function obeys dynamical scal- 
ing, i.e. (y{h)y(t 2 )) = tf*(ti/* 2 )- Let T = hit and 

x{T) = y(t)/ (y 2 (t)) . Then x(t) is a Gaussian station- 
ary process with zero mean, i.e. its correlation function 
is translationally invariant, {x{T\)x{T2)) — A(T2 — T%). 
Notice that A(0) — 1 by construction, which convention 
we shall adopt throughout this Communication (in con- 
trast to that of ref. H). If the persistence probability of y 
decays algebraically in t, then the persistence probability 
of x(T) decays as ~ exp(— 8T) for T — ► 00. 

The persistence probability may be expressed as the 
ratio of two path integrals, as follows M : 

f ^ n Dx(T)exp(~S) 
where 



S= - 



If dT x f dT2x(Ti)G(Ti,T 2 )x(T a ). (2) 
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Here G{T\,T2) is the matrix inverse of the correlation 
matrix (x(Ti)x(T 2 )) = A(T 2 - li). Notice that G is not 



simply a function of T 2 — T\ (unless we impose periodic 
boundary conditions). 

In MS this path-integral formalism was used to map 
the Markov process onto a quantum harmonic oscilla- 
tor in imaginary time, developing the perturbation the- 
ory in the formalism of quantum mechanics. We shall 
merely use path integrals as a convenient notation, per- 
forming all our calculations within the natural framework 
of stochastic processes. 

Let x°(T) be a stationary Gaussian Markov process, 
i.e. one defined by 



dx_° 
dT 



— = -vlx° + S(T), 



(3) 



where £ is a Gaussian white noise, with (£(T)£(T")) = 
2fj,d(T — T"). The noise strength has been chosen so that 
the autocorrelation function is A°(T) = cxp(— fiT). 

Suppose the process x(T) is perturbatively close to a 
Markov process, in the sense that G = G° + eg. Then 
we can expand the exponentials in the path integrals in 
(flh and re-exponentiate, so that to O(e) the numerator 
becomes 

J Dx(T)e~ s = J Dx{T) exp {-S° - | 7 (T) + 0(e 2 )) , 

(4) 

where the subscript c represents the constraint x(T') > 
(0 < T" < T) on the paths in the integral in the numer- 
ator of (|l|), and 



l{T) 



dTi / dT 2 g(T 1 ,T 2 )A° c (T 1 ,T 2 ), 
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where 



A° c (Ti,T 2 ) 



/ c Da;(T)a:(ri)a;(r 2 )f 
f c Dx(T)e-*> 



(5) 



(6) 



is the correlation function for the Markov process, av- 
eraged (and normalized) only over the paths consistent 
with the constraint C. The denominator in (pi) is given 
by an identical expression, except that A® is replaced by 
A , the unconstrained correlation function. 

By virtue of the constraint, AQ will not be strictly 
translationally invariant for finite T. In the limit T — > 
oo, however, the double time- integral in (|4|) reduces 
to T times an infinite integral over the relative time 
T 2 — T\, with A®(Ti, T 2 ) replaced by its stationary limit 
A^q{T 2 — T\). Similarly, g will be translationally invariant 
in this regime, giving 



7(T) 



(du/2Tr)g(Lo)A c (u;) , 



(7) 



where we have used the translational invari- 
ance to write the final result in Fourier 
space |L4J. Note that the zeroth-order result 



J x>0 Dx(T)eM~S )/ J Dx(T)exp(-S a ) is just the 
persistence probability of the stationary Gaussian 
Markov process x°(T), which decays as exp(— fj,T) as 
T — > oo. 

Using (|l|), (H) and (0), we find that the persistence 
exponent may be written in the form 

9 = lim -i log [P(x(T') > 0; < T' < T)] 

1 ^OO 1 

= M + e/ ^5(w) igM-i°M +0(e 2 ). (8) 



2tt- 

where the term in A°(ui) is the 0(e) contribution from the 
denominator in dlh, and we have exploited the u> — > — u> 
symmetry of the integrand. 

We now calculate Aq(T). The conditional probability 
Q(x,T\xo, 0) for the stationary Markov process may be 
obtained directly from (0): 

Q(x,T|xo,0) = 



2^(l-e-^ T ) 



1/2 



exp 



x — XQe 



-fiT\ 



2 ( 1 _ e -2pT) 



(9) 



The conditional probability Q + (x 2l T 2 \xi,Ti) that the 
process goes to (x 2 ,T 2 ), given that it started from 
(xi,Ti), without x ever being negative is given by the 
method of images: 

Q + (2|l) = Q(x 2 ,T 2 \xi,T 1 ) - Q{x 2 ,T 2 \-x 1 ,T 1 ), (10) 

where we have adopted an obvious shorthand notation 
for the arguments of Q + . 

To calculate the joint probability P + (xi,Ti;x 2 ,T 2 ) 
that the process passes through X\ at T\ and x 2 at T 2 , 
averaged only over paths where x(T) is always positive, 
we consider a path starting at (xi,Ti) and finishing at 
(Xf,Tf), passing through (xi,Ti) and (x 2 ,T 2 ) without 
ever crossing the origin. Then the required stationary 
limit is 



P + (x ll T 1 ;x 2 ,T 2 ) 



Ti- 



hm 9llf^l. (11) 

-oo.T^oo Q+{f\i) K ' 



The Markov property means that we can decompose 
t(/;2;JJt) = Q+(/|2)Q+(2|l)Q+(l|z). Using (§) and 
we find 



Q+(/;2;llt) 

P + (x 1: 0;x 2 ,T) 



(l-e-^ T y 1/2 x lX2 e» T x 



exp 



2(1- e-^ T ) 



sinh 



X\X 2 



2 sinh /iT 



(12) 
(13) 



It is now straightforward to evaluate the autocorrelation 
function: 

A° C {T)= dxi / dx 2 x 1 x 2 P + {x 1 ,0;x 2 ,T) (14) 
Jo Jo 



3 (1 - e- 2 » T ) 1/2 



+ (e 



M T 



2e"" T ) sin" 1 e-» T 



(15) 



Eq. (||) for 9 can now be expressed as a real-time in- 
tegral as follows. We first write A(T) = A°(T) + ea(T), 

r - i _1 

and we note that in Fourier space A(ui) = G(lu) — 

G°(w) + eg{w). Using A" = exp(-^T) gives g(w) = 
—a{uj){uj 2 + /i 2 ) 2 /4/i 2 . Inserting this in (H), and trans- 
forming to real time, gives 

/>oo 

e = ^-J^J o dTa(T)(n 2 ~d 2 T ) 2 [A°(T)-A°(T)] 



/i 



1 



2fi 



a(T) 



[1 - exp(-2/iT)] 



3/2 



dT 



(16) 



The final result is remarkably compact. Since ea(T) is 
just the perturbation to the Markov correlator A°(T) = 
e _MT , the normalization A(T) = 1 forces a(0) = 0. This 
is sufficient to converge the integral in (|l q ) p rovided a(T) 
vanishes more rapidly than T 1 / 2 . Eq. (|16[) has recently 
been used to calculate persistence exponents for inter- 
face growth in a class of generalized Edwards- Wilkinson 
models Q. 

As was remarked earlier, the problem of non- 
equilibrium critical dynamics is Markovian to first or- 
der in e = 4 — d. In the thermodynamic limit the 
global order parameter is Gaussian because, at time t, 
it is the sum of [L/£(t)] d (essentially) statistically in- 
dependent contributions, where L is the system size 
and £ ~ t 1 ^ is the length scale over which criti- 
cal correlations have been established. Corrections to 
the Gaussian distribution can be expressed in terms of 
higher cumulants of the normalized total magnetization 
M(i)/(M 2 (i)) 1 / 2 . Using the translational invariance of 
the system with respect to space it is easy to show that 
for large L the 27V-point cumulant is smaller by a fac- 
tor (t 1 ' z /Ly N ~ l ' d compared to the Gaussian part of 
the 2./V-point correlation function. The perturbative ap- 
proach discussed in the first part of this Communication 
can therefore be applied. To calculate the lowest non- 
Markovian term in 9, we need to calculate the autocor- 
relation function of the total magnetization M(t) to or- 
der e 2 , i.e. we need to calculate the autocorrelation func- 
tion A(t u t 2 ) = (Af(i 1 )Af(i 2 ))/(M 2 (i 1 )) 1 /2( M 2 (i2)) i/2 5 

which in the scaling regime depends only on the ra- 
tio t 2 /ti. The necessary techniques of dynamical field- 
theory, incorporating the extra renormalization associ- 
ated with the random initial condition (and responsible 
for the nonequilibrium exponent A) , have been developed 
by Janssen et al |15|,|l6|]. We first consider 'model A' dy- 
namics p7| for a nonconserved 0{n) model (where 'per- 
sistence' is associated with a given component of the or- 
der parameter for n > 1). The calculation is straight- 
forward in principle (but algebraically tedious), and the 
final result is (with T = ln(i 2 /£i)) 



A{T) = 



- P -^T 



1 



i^f^ T > + °e 3 > 



, (17) 



where fj, = (X — d+1 — r//2)/z from the one-loop calcu- 
lation 1 11 1 (equivalent to the 'Markov scaling relation'), 



and 



F A (x) 



- In - [2 ln(2a;) + (x - 1) ln(a; - 1) - (x + 1) m(o> + 1)] 



-2 (In 2) - 

+(x + l)ln 
-(x+ l)ln 
-(3a; 4-1) In 

2 
-(s-l)Li 2 

-2(s + l)Li 2 

+ (x + l)Li 2 



71 

x - 



-41n2- (x-l)ln 



x ■ 



1 



2x 

x+l 
2x 
3x 4 



+ (a; - 1) In 

, ,3x4-1 

In 



2x 
x-l 

2x 



In 



3a; -1 

2a; 



2x 



1 



2.7' 

( 3x - 1 
V 2x 
x- 1 

2a; 

X + l 



+ (3a; - 1) In 
, (s + 1) 



3a;- 1 
2a; 

3a; 



In 



2.C 



4.7' 

2a; 
3x4-1 



+ (x+ l)Li 2 
- 2(x - l)Li 2 
-(x-l)Li 2 



1 



2.7: 
x - 



1 



4a; 

2x 

3a; -1 



and Li 2 (a;) = — L dt ln(l — t)/t is the dilogarithm func- 
tion. The function FA(e T ) is a bounded, monotonically 
increasing, function of T in (0, 00). It vanishes as TlnT 
for T — > [satisfying the requirement for convergence at 
T = of the integral in @], while F(oo) = 0.057622 . . . 
The non-Markov nature of the process M(t) at order e 2 
follow from the fact that, at this order, A(T) is no longer 
a simple exponential. Substituting a(T) = A(T) — e~ MT 
from @ into Eqn. ©, using fi = (1/2) + 0(e), we find 



(after some algebra) 



3(n 4- 2) 



e 2 a 



(18) 



, 4(n4-8) 2 
where 
a = 4 (\/2 - 2V3 4- VE) 4-8\/2hi2-4fV2-l)ln3 
-2 (l + 2V2) In (3 4- 2V2) - 14 In (5 4- 2\/6 

+ 101n(7 + 4V3)+8^1n ^ + ^-^ ) 

-4V21n f 2 f- 2 + fl 
\2V3-2-V2J 

= 0.271577604975 . . . 

This result can be compared with recent simulation 
data for the Ising model in 2 (ll]jl|@ and 3 Jl| di- 
mensions. For d = 2, using A = 1.585 ± 0.006 pl|, and 
rj = 1/4 (exact) gives \iz = 0.460 ± 0.006. Ignoring non- 
Markov corrections, one would obtain 9z — /iz, smaller 



than the measured value 9z = 0.505 ± 0.020 (the finite- 
size scaling method used in |11| naturally determines the 
combination 9z pQ]). The non-Markov correction fac- 
tor in (|J) is, for n = 1, (1 + 0.0075438 . . .e 2 ) ~ 1.030 
for e = 2. The 'improved' estimate for Qz becomes 
0.474 ± 0.006, closer to, but still somewhat smaller than, 
the numerical estimate. 

For d = 3, one has z = 2.032 ±0.004, A = 2.789 ±0.006 
@ and n = 0.032±0.003, giving fi = 0.380±0.003. Mul- 
tiplying by the non-Markov correction factor for e = 1, 
i.e. 1.0075, gives 9 = 0.383 ± 0.003, compared to the nu- 
merical result 9 ~ 0.41 Q. A direct expansion to order 
e 2 , using the known expansions for z, A, and rj, gives (spe- 
cializing to n = 1) 6 = 1/2 - e/12 + (a - 21n3)e 2 /72 - 
2e 2 /81 + <3(e 3 ), i.e. (9 ~ 0.365 for d = 3, slightly lower 
than that obtained using the best numerical estimates of 
z, A and r\ and only using the e-expansion for the non- 
Markov correction. 

A similar approach can be applied to 'model C critical 
dynamics fll7| , in which a nonconserved order parame- 
ter field is coupled to a conserved density. In this case, 
one obtains non-Markovian corrections already at order 
e. The autocorrelation function is given by (for n = 1) 



A(T) = exp(- M T) [l - - F c (e T ) + 0(e 2 )\ (19) 

x — 1 

Fc(x) = In 2 h a: In x 

2x 

^ln(^l)-^lln( I + l). (20) 



The work of SC and AB was supported by EPSRC 
grants GR/K53208 and GR/J24782 and the work of KO 
was supported by grant Oel99/l-l of the DFG. 



Again, F c (e 
F c (oo)=ln2-Jj2. 
from (|19|) into ( |16[) g: 



vanishes like TlnT for T 

Inserting a{T) = A{T) 

1S() into ( |16[) gives 



7-4^2 



0{e 2 ) 



-> 0, while 
exp(-^T) 



(21) 



where /i = (A — d + 1 — rj/2)/z as before, but now the 
dynamical exponents 2 and A take their model-C values 

In summary, we have computed to order e 2 the persis- 
tence exponent 9 for the global order parameter M(t) of 
models A and C. At this order, the dynamics of M(t) are 
non-Markovian, and 9 is a new exponent, not related to 
the usual static and dynamic exponents. The calculation 
was performed by expanding around a Markov process, 
using a simplified form of the perturbation theory intro- 
duced by Majumdar and Sire. 
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